Hydrino like states in graphene and Aharonov-Bohm field 
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We study the dynamics of fermions on graphene in presence of Coulomb impurities and Aharonov- 
Bohm field. Special emphasis is given to the formation of hydrino like states and its lifting of 
degeneracy due to the presence of AB field. The flux of the AB field can be tuned to make the low 
angular momentum hydrino states stable against decay. The zero limit physics of the two coupling 
constants etc and $ involved in the system is discussed. 
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Graphene, a 2-dimensional form of allotrope, made up 
of single layer of carbon atoms, has been a platform to 
study the various exotic features [l| of fermions which is 
usually not seen in most condensed matter systems. The 
exotic behavior is attributed to the unique lattice struc- 
ture of the allotrope, which can be understood from the 
fact that the different layers of graphite are so loosely 
bound that it slides one above another. The bonding be- 
tween the atoms of the same layer are strong enough to 
form a stable structure. It can be seen when struck by 
a pencil lead on a paper. The mark on the paper can 
be considered as a thin layer of carbon atoms. The re- 
cent experimental formation of graphene in 2004 0] has 
made it possible to test several features in this system. 
For example, in semiconductor or in other material, the 
energy-momentum relation for the electrons with effec- 
tive mass m* are given by E = p 2 /2m* . But for mono- 
layer graphene without any external perturbation or dis- 
tortion, the dispersion relation for the electrons near a 
Dirac point are given by E — ±cq\p\ with velocity cq. 
This linear dispersion relation was a clue that the charge 
carriers in graphene could be described by massless Dirac 
equation. However, mass term for the fermions could 
arise if somehow the two sub lattices of the graphene are 
made distinguishable. 

One of the exotic behaviors is that, the mobility of 
charge carriers in graphene is approximately indepen- 
dent of the carrier density. The constancy of the mo- 
bility means that the conductivity in graphene is pro- 
portional to the carrier density. The list of unusual phe- 
nomena in graphene is even extended in quantum Hall 
effect, which can not be explained by non-relativistic 
quantum mechanics. The integer quantum Hall effect 
occur at half integral values of the filling factor, which is 
therefore termed as anomalous [l], 0] . This effect is com- 
mon in massless monolayer graphene and massive bilayer 
graphene. The Landau levels which are important to de- 
scribe quantum Hall effect are also different from what 
is seen in non-relativistic quantum mechanics. In con- 
stant magnetic field B background perpendicular to the 
plane, the electron motion in metal or semiconductor is 
quantized as E — hto*(n + 1/2), n = 0, 1,2, ... The holes 
also have same behavior but with negative energy. Here 



H is the reduced Plank constant and w* = eB/m* is the 
cyclotron frequency. But graphene has completely dif- 
ferent Landau levels In monolayer graphene, it is 
given by Eq = sgn(n)cc y/2eh\n\B . The occurrence of 
the Eg = states are the consequence of chirality 0| 
of the Dirac fermions. The degeneracy of the Eg = 
state is however same as the Eg ^ state. On the other 
hand, for bilayer graphene the spectrum is of the form 
E G = sga(n)Hu Gy /\n\(\n\ + l),n = 0,±1,±2,.... The 
Eg = levels are still present but this time the degener- 
acy is double than that of Eg ^ levels. 

All the above features are indicative of the fact that the 
dynamics of charge carriers in graphene is described by 
Dirac fermions. The tight binding model || of graphene 
indeed show that the dynamics of electron and hole are 
dictated by (2+l)-dimensional Dirac fermion in the con- 
tinuum limit of the lattice spacing. Graphene is thus 
considered to be a condense matter analog of QED (see 
the introduction of Q). However, it has lot of differences 
with the standard QED. Depending upon the situation, 
fermions could have mass or zero rest mass as mentioned 
before. The velocity of the zero rest mass fermions are 
C G ~ c/300, which is much less than the light velocity 
c in vacuum. Due to this small Fermi velocity cg, the 
fine structure constant is much higher ag « 300a. Ob- 
viously the perturbation scheme of QED seems to be not 
applicable to graphene system due such a high coupling 
constant olg- 

In this work, we report another possibility of exotic 
phenomena of massive fermions on graphene, which is 
the formation of hydrino like states. But due to the 
high coupling constant olg, the so called hydrino states 
may have complex eigenvalues for low angular momen- 
tum values, which will make it dissipative. We therefore 
consider a Aharonov-Bohm (AB) field, which will rescue 
the low angular momentum states and will form stable 
bound states. Hydrino state is the tightly bound states 
of Hydrogen atom with binding energy Eb > 13.6eV, 
where normal ground state energy for Hydrogen atom is 
Eg = —13.6eV. Experimental study has also been re- 
ported recently in support of this hydrino states (see 
and its refs.). Theoretical model of hydrino states can 
be found by solving Klein- Gordon equation in 3 space 
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dimensions. Recently it has been shown that, in (2 + 1)- 
dimensions, Dirac equation may give rise to hydrino like 
states [8|. Besides the hydrino states, the large fine struc- 
ture constant qq of graphene also may make the normal 
bound state eigenvalue complex. Our inclusion of a AB 
field will also make the normal bound states stable. 

The existence of deeply bound states can be explained 
in terms of a consistent mathematical formulation known 
as self-adjoint extensions [9(. In this method the com- 
monly used Dirichlet boundary condition at origin of the 
wave- function, that the wave-function is zero at origin, 
is generically extended to a less restricted boundary con- 
dition. Usually the wave-function at the origin is con- 
sidered to have a certain power law behavior without 
violating the square-integrability condition. This larger 
domain then allows other states to immerse out as a vi- 
able solutions. 

We start our discussion with the fermions on graphene 
in presence of a gauge potential A^, described by a two 
component Dirac equation 

(i$ + A-M)i>(t;r,cl>) = (1) 

We set the unit in which Fermi velocity cq « 10 6 m/,s = 1 
and H = 1. Here <j? = 7 M 9 M , A = "f^A^ and we take the 
representation of the gamma matrices in terms of Pauli 
matrices as 7 = a 3 , 7 1 = ia 1 and -f 2 = ia 2 . Since we 
are considering Coulomb impurities and Aharonov-Bohm 
field, our gauge potential A^ is the Coulomb potential 
Aq = ac/r and AB field A$ = -$/r. The mass M of 
the fermion could generate due to the external perturba- 
tion, which makes the two sub-lattices of the honeycomb 
lattice structure of graphene distinguishable or by other 
means [Io|. The the hopping of the electrons between 
two sub-lattices makes the energy difference nonzero. In 
order to solve ([1]) we need to know the condition under 
which the radial hamiltonian extracted from (flj is sym- 
metric. This condition can be easily found to be 

lim r<j)\{r)i(j 2 4> 2 {r) = 0, (2) 

where cf>i and $2 are two component radial spinors. To 
construct a domain for our problem we need to construct 
it in such a way that it respect the constraint ([2]) and 
at the same time it becomes equal to its adjoint domain. 
We mentioned in the introduction that, the mathemat- 
ical formulation which serve this purpose is known as 
self-adjoint extensions and it has been recently discussed 
in case of graphene with Coulomb impurities in and 
a generalized boundary condition for the bound state so- 
lution has been obtained. 

We do not go into detail analysis of it, rather we give 
here an intuitive argument to show that the boundary 
condition can be parameterized by U(l) group. The short 
distance behavior of the two independent radial wave- 



functions %px (r) and ip2 (f) has the form 
lim^xfr) - Ar- 1/2+i 

r—>0 

lim-i/^r) ~ Br- 1 / 2 ^, (3) 

where A and B are the two component constants, £ = 
v / (m + $) 2 -a G , and m = ±1/2, ±3/2, .... For $ = 0, 
<j3j> is reduced to the short distance behavior of radial 
function of [8(. Note that both the solutions are square- 
integrable at r = if 

< 2\Re(0\ < 1 (4) 

Technically, the point r — > is called quantum mechan- 
ically incomplete Q and the potential corresponding to 
the radial Hamiltonian is called in limit circle case. On 
the other hand only one solution is square- integrable near 
r — > 00. So the potential is in limit point case at infin- 
ity. This means, there can have 1-parameter family of 
self-adjoint extensions. Since both solutions are square- 
integrable at r — > 0, then any linear combinations of the 
two will also be square-integrable. We take a generalized 
boundary condition of the form 

lim Mr) = Cx(e)r- 1/2+e + C 2 (0)r~ 1/2 ~ S , (5) 

r—>0 

where the constants, C\(9) and €2(6), which are now de- 
pendent on the free parameter 9 of the U(l) group, be 
such that ([5]) satisfies the constraint ([2]) . Then 9 is called 
the self-adjoint extension parameter. The normal bound 
state solutions can be obtained by setting the coefficient 
of the singular term of (O to zero, i.e., 02(d) = 0. This 
gives the bound state solution of the form 

E n>m = Msga(a G )/]l 1 + a 2 G /(n +^[(m + $f- a 2 G f.(6) 

Note that for $ = the above solution reduces to the 
result d, [l3j]. Now comes the exotic solutions of the 
problem, which can be obtained by setting the 1st term 
of ([5]) to zero but keeping the singular term. In this case 
the solutions are square-integrable and the bound state 
eigen- value is given by 

E n . rn = Msga(a G )/^ 1 + a 2 G /(n - yf (m + $) 2 - a 2 G ) 2 ,(7) 

which reduces to the result of [8] for $ = 0. This is 
the so called hydrino like states. Note that in presence 
of the magnetic flux, the angular momentum degeneracy 
is lifted. As is evident the states with spectrum ((7|) are 
tightly bound than the states with eigenvalues ([6]). we 
have two parameters in our system, one is fine structure 
constant ag and other one is the flux whose limiting 
behavior can be investigated. We see that the limit (f> — » 
0^ is smooth and from ^ and (0 we get 

lim E n m = lim E n m = lim E n , m (8) 
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But the similar limit for ag shows that 



lim E„ 



7^ lim 



En r, 



^ lim E n r . 

a G =0 



(9) 



Note that the first inequality of (J9j) can be obtained from 
© or (|7|), but to get the second inequality we need to 



compare our result with 12 1 



To summarize, we considered the quantum dynamics 
of fermion on graphene in presence of Coulomb impurities 
and AB field background. Possible generalized boundary 
condition has been established based on intuited argu- 
ments. The fermions are trapped in the combined poten- 
tial. We showed that the deepest bound state solutions, 
known as hydrino like states, are nothing but a conse- 
quence of a specific choice of boundary condition. It can 
be rigorously shown by self-adjoint extensions. The limit 
$ — > is shown to be smooth but the limit of the fine 
structure constant is not smooth so can not be used to 
get the dynamics of the fermions in only AB field back- 
ground. Due to large value of the fine structure constant 
qg of graphene, the low angular momentum states E n<m 
may become complex when only Coulomb field is present. 



But the presence of the AB field may rescue it and make 
the states stable. 
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